Distribution of Cold dust in Orion A and B 
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ABSTRACT 

Large scale far-infrared (FIR) observations of the Orion complex at 205 and 138 /mi are presented 
with an aim of studying the distribution of cold (<25 K) dust. The maps in these FIR bands extend over 
-3600 sq. arcmin and cover regions around OMC-1, 2, 3 in Orion A and NGC 2023 and NGC 2024 in 
Orion B. Some limited regions have also been mapped at 57 ^m. A total of 15 sources in Orion A and 14 
in Orion B (south) have been identified from our FIR maps. Dust temperature distribution in both Orion 
A and Orion B (south) have been determined reliably using the maps at 205 & 138 /im obtained from 
simultaneous observations using almost identical beams (1.'6 dia). These temperatures have been used to 
generate map of riso, the optical depth at 150 /im, for the Orion B region. The coldest source detected 
is in OMC-3 and has a temperature of — 15 K. The diffuse FIR emission in the different sub-regions is 
found to vary between 25 % to 50 % of the total FIR emission from that sub-region. 

Subject headings: infrared: ISM: continuum - ISM: individual (OMC-1, OMC-2, OMC-3, NGC 2024, NGC 2023) 



1. Introduction 

In this paper we present a large scale high angular 
resolution (~1') study of the continuum emission at 
205 and 138 /im from the two Giant Molecular Clouds 
Orion A and Orion B (south) in the Orion complex. 
The mapped region of Orion A includes OMC-1, 2, 
3 and IRAS 05327-0457 (NGC 1977). In Orion B 
we have mapped mainly the southern portion covering 
NGC 2023 and NGC 2024. Limited regions including, 
OMC-1 and 2 in Orion A and NGC 2024 in Orion B 
have also been observed at 57 /im. 

The southern part of the Orion cloud complex, also 
known as Orion A, contains star forming cores which 
are distributed along the north-south direction (integral 
shaped ridge) and named OMC-1, 2 and 3 from south 
to north. Observations at sub-millimeter/far-infrared 
wavelengths have shown that both warm and cold dust 
with quite different spatial distribution contribute to 
the emission from the Orion molecular clouds (Thron- 
son et al. 1986; Chini et al. 1997). Cold dust, with 
temperature — 20-25 K, accounts for bulk of the mass 
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of dust in the Orion clouds and emits significantly only 
at wavelengths longer than 100 /im. Earlier observa- 
tions of Orion A concentrated on the emission from 
the hot dust viz., mapping of OMC-1 at 50-100 /im 
by Werner et al. (1976) and of M43 at 37-160 /im 
by Smith, Harper, & Loewenstein (1987). Later ob- 
servations have mostly explored the emission from the 
cold dust in this region. The most recent of these ob- 
servations was by Johnstone & Bally (1999) covering 
OMC-1, 2 and 3 and extending over an area of 50' x 10' 
with an angular resolution of — 7"5 at 450 /im and 14" 
at 850 /im. Lis et al. (1998) had also mapped a region 
of 6'x 18' covering OMC-2 and 3 at 350 /im with an 
angular resolution of 11". However, the pioneering 
work in detecting a number of cold sources in OMC-2 
and an extremely cold one (~ 10 K) in OMC-3 was by 
Chini et al. (1997), who mapped a region of — 15'x5' 
with an angular resolution of — 1 1" at 1300 /im. Most 
of these condensations are considered to be protostel- 
lar candidates because of the associated outflow activi- 
ties (Castets & Langer 1995). Observations in molec- 
ular lines like C ls O by (Castets & Langer 1995, ;and 
references therein) have revealed that the gas in OMC- 
2 and OMC-3 have kinetic temperatures equal to 24 
K and 19 K respectively, both much smaller than that 
in OMC-1 (70 K). These observations also identified 
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both, OMC-2 and OMC-3, to be regions of high col- 
umn density. C ls O observations have also detected 
several moderate to strongly collimated bipolar out- 
flows in OMC-3, which suggests recent star formation 
activity. The very low temperatures of sources in both 
OMC-2 and 3 along with the evidence for star forma- 
tion, indicate very early stages of evolution of, prob- 
ably intermediate mass stars. In this paper we have 
explored the distribution of dust emission and temper- 
ature of the Orion A region in a wavelength regime 
(205 & 138 /im in particular) where no high angular 
resolution measurements exist in literature. 

The Orion B cloud contains 4 major sites of star for- 
mation viz., NGC 2071 & NGC 2068 in the north and 
NGC 2023 & NGC 2024 in the south. We have ob- 
served only the regions around NGC 2024 and 2023. 
The brightest emission source in the region is NGC 
2024, which is a region of active star formation. The 
1300 /im observations over a region of <~ 5'x5' (with 
an angular resolution ~ 30") by Mezger et al. (1988) 
have revealed six condensations in NGC 2024, aligned 
in the north-south direction. These condensations also 
have masers and HII regions associated with them. 
The nature of these condensations was an issue of de- 
bate which was resolved by Mezger et al. (1992). 
From their observations at 870 & 1300 /im they con- 
cluded that these condensations are indeed isothermal 
protostars. Lada, Bally, & Stark (1991a) presented 
a survey of Orion B (Lynds 1630) molecular cloud 
in CS(2-1) line ( 1 '8 angular resolution) and identified 
the regions including NGC 2023 and NGC 2024 to be 
bright in CS, thus indicating high density of the region. 
In a follow-up program Lada et al. (1991b) surveyed 
Orion B in the K (2.2 /im) band and found that more 
than 50% of the NIR sources are associated with the 4 
identified star forming regions. In a recent high an- 
gular resolution (22-30") survey of Orion B (south) 
at 1300 /im extending over ^80'x68', Launhardt et 
al. (1996) have identified some of the dust emission 
peaks around NGC 2023 and NGC 2024 with the CS 
cores detected by Lada, Bally, & Stark (1991a). Of 
the existing FIR observations of NGC 2023 and NGC 
2024, the IRAS observations in the first place trace 
relatively warmer dust (> 30 K) and secondly have 
angular resolutions varying with wavelength, so that 
calculating the temperature distribution from the ob- 
served flux densities is not straightforward. On the 
other hand, the 1300 /im observations, though of high 
angular resolution, trace extremely cold dust. This has 
motivated us to carry out simultaneous high angular 



resolution far-infrared observations at 205 and 138 /im 
with identical beams so that the density and tempera- 
ture distribution of dust emitting significantly in this 
intermediate regime can be reliably quantified. 

The rest of the paper is organized as follows: Sec- 
tion 2 describes the observations and Section 3 con- 
tains the results and discussion of the same. 

2. Observations 

2.1. The 205 & 138 /im maps 

Orion molecular clouds A and B (south) were ob- 
served using a two-band far-infrared photometer sys- 
tem (PHT-12) at the Cassegrain focus of the TIFR 1- 
meter (f/8) balloon-borne telescope. The far-infrared 
telescope was flown from the TIFR Balloon Facility, 
Hyderabad, in central India (Latitude = 17.°47 N , Lon- 
gitude = 78P57 E ) in November, 1995. Details of the 
telescope and the observational procedure have been 
given by Ghosh et al. (1988). The photometer (PHT- 
12) consists of 12 composite silicon bolometers, each 
having a FOV of 1.'6 and arranged in a 3x2 array for 
each band. These arrays are cooled to 0.3 K by liquid 
3 He. The sky was chopped at 10 Hz by wobbling the 
secondary mirror. The chopper throw was 4f2 along 
the cross-elevation axis. The same region of the sky 
was viewed simultaneously in two bands. The effec- 
tive wavelengths for the two bands for a A" 2 emissivity 
law and a temperature of 35 K are 205 /im and 138 
/im. Saturn was observed for absolute flux calibration 
(Ghosh et al. 1988) as well as determination of the in- 
strumental Point Spread Function (PSF) including the 
effect of sky chopping. 

The simultaneous mapping in the two FIR bands 
was carried out by using raster scans in which the 
telescope was continuously moved along the cross- 
elevation axis and stepped in elevation. The chopped 
FIR signals were gridded in a two dimensional sky ma- 
trix (elevation x cross-elevation) with a cell size of 
0'3 x 0'3 . The observed chopped signal matrix was 
deconvolved using an indigenously developed proce- 
dure based on the Maximum Entropy Method (MEM) 
similar to that of Gull & Daniell (1978) (see Ghosh 
et al. 1988, for details). The FWHM sizes of the de- 
convolved maps of the point-like source (Saturn) are 
found to be 1 .'6 x 1 '.9 and 1 .'6 x 1 .'8 in the 1 38 and 205 
/im bands respectively. An optical photometer was 
used at the Cassegrain focal plane to improve the ab- 
solute positioning of the telescope during the flight. 
In addition for all the observations presented here, 
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the signal from this optical photometer measured dur- 
ing the far-infrared scans have been used during of- 
fline analysis to further improve the aspect informa- 
tion. The absolute positional accuracy achieved in the 
FIR maps is ~ 0.'5 . 

2.2. The 57 /im map 

Limited regions of both Orion A (OMC-1 & 2) 
and Orion B (NGC 2024) were also observed at 57 
/im in an earlier flight (November, 1990) of the TIFR 
balloon-borne telescope, using a different photome- 
ter (PHT-2). This two band photometer, PHT-2, was 
based on two liquid 4 He cooled Ge(Ga) bolometers. 
The FIR passbands were 45-70 /jm and 110-160 /mi 
with effective wavelengths 57 /mi and 150 /im for an 
input source spectrum defined by a greybody of tem- 
perature 35 K and A" 2 emissivity. A 2'A (dia) re- 
gion of the sky was simultaneously viewed in both 
the FIR bands. More details of PHT-2 can be found 
in Verma et al. (1994). Here, only the observations 
at 57 /im are presented, since the other band is very 
similar to the shorter waveband of PHT-12 which has 
better sensitivity. The PHT-2 was calibrated by ob- 
serving Jupiter. Some preliminary results from obser- 
vations of the Orion region have been presented by 
Das et al. (1993). The final results are presented 
here. The achieved angular resolution in the image 
processed maps is <~1?5 as determined from the ob- 
servations of Jupiter. 

2.3. The HIRES processed IRAS maps 

To supplement our balloon-borne observations, we 
have used the IRAS survey data in all wavebands (12, 
25, 60 and 100 /im) for the regions mapped by us in 
Orion B (south). These data were HIRES processed 
(Aumann, Fowler, & Melnyk, 1990) at the Infrared 
Processing and Analysis Center (IPAC 1 , Caltech) for 
improving the angular resolutions of the raw maps. 

3. Results & Discussion 

3.1. Analysis of FIR data 

Figure 1 and Figure 2 respectively show the regions 
(indicated by the dotted boundary) in Orion B (south) 
and Orion A mapped at 205 and 138 /mi. Regions ex- 
tending over ~ 42'x38' in Orion A and ~ 38'x53' 

'IPAC is funded by NASA as part of the IRAS extended mission 
program under contract to JPL. 



in Orion B (south) have been observed at both these 
FIR wavelengths. The 57 /tm maps for both sources 
cover relatively smaller regions. Reliable dust temper- 
ature T( 138/205) maps for selected regions have been 
generated utilising the simultaneity of observation and 
nearly identical beams at 205 and 138 /on. The dust 
emissivity law has been assumed to be e„ oc v& and the 
method used is briefly described in Appendix 1. For 
most star forming regions the emissivity exponent at 
longer wavelengths (> 100 /im) is found to be between 
1 .5 to 2.5 and the calculated properties of most com- 
mon types of dust grains (Draine & Lee 1984) show 
the same at these wavelengths. Based on these facts we 
have chosen (3 = 2 for deriving the T( 138/205) distri- 
bution. Considering the uncertainties in the measured 
flux densities and the wavelength regime probed, the 
temperature computed has an uncertainty of ~ ±2K 
between 14 K and 50 K and of ~ ±5K between 50 and 
70 K. The intensity ratio of our FIR bands is rather in- 
sensitive to temperatures below 10 K and above 70 K. 
Above 70 K the estimated temperature has an uncer- 
tainty of ~ ±8K. The morphology of the spatial distri- 
bution of Td is insensitive to the assumed value of (3. 
Only the numerical values of T<j depend on (3, e.g., a 
temperature of 15 K for (3 = 2 corresponds to 17 K for 
/3 = 1 and Td = 30 K for (3 = 2 corresponds to 50 K for 
[3=1. Distribution of ti 5 o, the optical depth at 150 /im 
for Orion B has also been generated using these dust 
temperatures and (3 = 2. 

We have adopted a well-defined extraction algo- 
rithm to identify discrete sources in these maps. This 
algorithm is based on identification of a local maxi- 
mum followed by a critical study of the growth curve 
of the integrated flux densities in consecutive annular 
rings with the maximum at the center. The reliabil- 
ity of the sources detected in each map has been en- 
sured by the following procedure. The entire observa- 
tion for each region has been split into two indepen- 
dent datasets (by considering only the alternate scan 
lines in each raster) and processed (deconvolved) sep- 
arately to generate two maps. The peaks detected in 
both these maps within a predetermined positional er- 
ror have been used to quantify the dynamic range and 
also to validate the detected sources. The final maps 
presented combine all data and positions of the "con- 
firmed sources" in the composite maps have been de- 
termined. These sources have signals more than 10 
times the estimated noise level of the respective maps. 
In this paper the method of presenting the sources is 
as follows : Peaks identified in any two TIFR wave- 
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bands are associated if they lie within 1'. For associ- 
ation with IRAS (HIRES) sources the same parameter 
has been taken to be 1.'5. The coordinates quoted for a 
particular source are for the position of the peak at the 
longest (TIFR) wavelength at which it is detected. 

3.2. Orion A 

Figure 2 shows the 205 and 138 /im intensity dis- 
tribution of the mapped region in Orion A. The maps 
include the Orion Molecular Clouds 1, 2, 3 and the 
source IRAS 05327-0457. Figure 3 shows the inten- 
sity map for a limited region in Orion A at 57 /im. Ta- 
ble 1 presents the coordinates and flux densities (in a 
circle of 3' dia) for all the discrete sources detected in 
the mapped regions of Orion A at 205 and 138 /mi. In 
the same table we also present the longer wavelength 
associations of these sources along with the measured 
flux densities at 1300 /im (Chini et al. 1997). For 
the discrete sources detected at 138, 205 and 1300 /im 
the dust temperature (Td) and the emissivity exponent 
(/?) have been determined by fitting a single tempera- 
ture greybody spectrum. Td and (3 for A8 have been 
determined in a similar manner by using the flux den- 
sities measured at 450, 790 and 1100 /zm by Gold- 
smith, Bergin, & Lis (1997) along with those from 
the present work. For all other sources with flux den- 
sities measured only at 205 and 138 /im Td has been 
determined assuming (3 = 2. For a source not detected 
in any one of the wavebands (marked by f in Table 1), 
we give the flux density obtained by integrating at the 
position of the source as detected in the other wave- 
band^). Table 2 presents the positions, flux densities 
and associations of these sources detected at 57 /im. 
The source IRAS 05327-0457 has been studied in de- 
tail in the far-infrared by Mookerjea et al. (2000); here 
we concentrate mainly on OMC-1, 2 and 3. 

Figure 4 shows the temperature distribution as cal- 
culated from the observations at 205 and 138 /im for 
OMC-1 and 2. The temperature in OMC-1 is as high 
as <~ 70 K close to the Orion nebula. The tempera- 
ture of the OMC-1 and the bar region as calculated by 
Werner et al. (1976) from the ratio of 50 and 100 /im 
flux densities is found to be between 55 and 85 K. For 
the same region, the temperature calculated from our 
205 and 138 /im fluxes is between 45 and 70K. This is 
consistent with the fact that the longer wavelength ob- 
servations presented here are tracing colder dust com- 
ponent as compared to those traced by the 50 and 100 
/im observations. The temperature drops to an average 
of - 25 K in OMC-2 and is as low as 15 K in OMC-3 



(not shown in Figure). The temperature map has been 
generated by assuming the dust to be optically thin but 
explicit calculation of optical depth (at 150 /im) very 
close to the global intensity peak of OMC-1, shows the 
regions to be optically thick. In such cases, the temper- 
atures quoted would be underestimates by as much as 
5K. However, for OMC-1, 2 and 3, the basic morphol- 
ogy of the isotherms remains unaltered. 

Orion A is an ensemble of cores with a variety of 
temperatures, physical sizes and masses. This variety 
in physical conditions makes it conducive to harbour- 
ing stars in different stages of formation and evolution 
and hence makes it astrophysically interesting. The 
far-infrared maps presented here also support this vari- 
ety by showing morphological characteristics strongly 
dependent on the dust temperature. While the 205 and 
138 /im maps tracing colder dust detect FIR emission 
from OMC-2 and OMC-3 very clearly, the 57 /im map 
fails to do so. 

OMC-1 (A8 in Table 1 and # 1 in Table 2) is the 
brightest source detected at all wavelengths. As men- 
tioned earlier, the submillimeter fluxes at 450, 790 and 
1 100 /im (Goldsmith, Bergin, & Lis 1997) have been 
used along with the 57, 138 and 205 /im fluxes to de- 
termine Tj and (3 for this source. The dust temperature 
obtained by this method is 68 K and the (3 is 0.5. If 
instead we do a similar fitting omitting the flux density 
at 57 /im we get Td =30 K and (3= 1.1. Thus the high 
temperature and low (3 as presented in Table 1 might be 
due to the effect of attempting to explain the emission 
over a wide range of wavelengths using a single tem- 
perature greybody spectrum. In reality there is prob- 
ably a cold diffuse component of dust as well which 
contributes to the emission at longer wavelengths. In 
the OMC-1 region the maps at 205 & 138 /im also de- 
tect the source A5. Although the position of this source 
at 205 and 138 /im differ by ^1'.5, we have associated 
the two, since there is no other nearby source of sim- 
ilar strength. We also note that the temperature of the 
source A5 is ~ 80 K, which is even higher than the 
temperature of A8. 

The peaks # 2 and # 3 detected at 57 /im (see Ta- 
ble 2) are not detected at the longer wavelengths. The 
source # 2 is associated with the HII region M43 and 
the intensity distribution around it is similar to that 
observed at 60 /im by Smith, Harper, & Loewenstein 
(1987). The source # 3 is due to the heating of the dust 
in the photodissociation region around the trapezium 
stars. Northern tip of the emission detected at 57 /im 
is clearly displaced to the east with respect to the inte- 
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gral shaped emission ridge detected in OMC-2 and 3. 
We can discern local enhancements in T( 138/205) map 
close to the positions of the peaks #2 and #3. These 
strongly indicate the existence of hotter dust towards 
the east. 

Of the observations presented here OMC-2 and 3 
are detected only at the longer wavelengths (205 and 
138 /im). In OMC-2 we detect two sources (A10 and 
All) at 138 /im and only one source (A12) at 205 /mi. 
As suggested by Table 1 the sources detected at our 
resolution are each associated with more than one sub- 
millimeter (1300 lira) sources. We note that whereas 
A10 and All are associated with two mutually exclu- 
sive groups of submillimeter sources, A12 is associ- 
ated with sources from both groups. The fact that we 
detect A10 and All only at 138 fmv and A12 only 
at 205 /im could be due to the presence of a diffuse 
cold component which reduces the contrast between 
the sources and the intermediate region. This could 
effectively smoothen out the smaller peaks (like A10 
and All) in this region , thus resulting in only one 
peak (A 12) at 205 /im. The flux densities at 1300 /urn 
of the detected sub-mm sources have been calculated 
using the peak fluxes and size of the individual sources 
(refer to Table 1 of Chini et al. 1997). For each FIR 
source the contribution at 1300 /im from all associ- 
ated sources are combined and presented in Table 1. 
For all sources detected in OMC-2 a single tempera- 
ture greybody spectral fitting gives 20K<Td<25K and 
1.6</3 < 2.0. Both these are consistent with the re- 
sults of Chini et al. (1997). Recent radio observations 
at 3.6 cm by Reipurth, Rodriguez, & Chini (1999) 
have detected sources coinciding with the peaks in dust 
emission. These associations are also presented in Ta- 
ble 1 . The radio emission further substantiates that the 
sources detected in OMC-2 though reasonably cold are 
not only regions of high density but also have YSOs 
embedded in them. 

We have detected one source (A9) in the OMC-3 
region. Positionally it is coincident with the sub-mm 
source MMS 6 (Chini et al. 1997), the latter being the 
brightest of all the sub-mm sources detected in OMC- 
3. The flux density at 1300 /im for the source A9 is 
estimated in the same manner as above. A greybody 
spectrum fitted to the measured flux densities gives Tj 
= 15 K and f3 = 2.0. A9 is thus the coldest source in 
Orion A, both Tj and (i estimated here are consistent 
with Chini et al. (1997) and Lis et al. (1998). 

The total far-infrared luminosity has been estimated 
using the temperature map (Fig. 4) for the region with 



OMC-1 and OMC-2 and the intensity maps of the 
same region. The dust emissivity exponent assumed 
for this calculation is 2.0. From the FIR observations 
presented here the total luminosity detected in OMC-1 
and 2 is ~3.5 x 10 5 L Q . For a similar region the lumi- 
nosity based on 60 and 100 /im observations as given 
by Stacey et al. (1993) is -5 x 10 5 L©. The diffuse 
emission at these wavelengths has been calculated by 
subtracting the contribution of the discrete sources and 
is about 25% of the total luminosity. 

The molecular cloud OMC-1 has long been the 
most extensively studied region in Orion A. OMC-2 
and 3 have been studied only in the recent past and 
have been found to be rich and interesting as a star 
forming region. Our observations provide useful con- 
straints on the intensity and temperature distribution of 
this region. 

3.3. Orion B (south) 

3.3.1. Region around NGC 2024 

Figure 5 presents intensity maps of the region 
around NGC 2024 at 205, 138, 100 and 57 /im. The 
map at 57 /im is for a smaller region, but brings out 
the main features of NGC 2024 at this wavelength. 
Table 3 presents the measured flux densities at the 
4 IRAS and 3 TIFR wavelengths for the sources de- 
tected in this region. The sub-mm and IRAS Point 
Sources associated with these FIR sources, along with 
the sub-mm flux densities measured at 1300 /im by 
Launhardt et al. (1996) are presented in the Table 3. 
In all 7 sources have been identified with sources in 
the IRAS Point Source Catalog (PSC) and 2 sources 
(B3 and B9) have been found to have sub-millimeter 
(1300 /tm) counterparts. For sources not detected in 
one or two of the TIFR bands (marked by f in Ta- 
ble 3) we have presented the flux densities obtained 
by integrating in a circle around the position of the 
source as detected at the other wavelength(s). We have 
fitted a single temperature greybody spectrum for all 
sources for which measurement of flux densities at 3 
or more wavelengths longer than 50 /im exist. The re- 
sulting Td and (3 are also presented in the same table. 
For the remaining sources Td is determined using (3 = 
2. We find that the emissivity exponent (/?) for some 
sources (B6, B10 and B12) are very small (< 0.5). 
In addition, for B6 and B12 the fitted Td is also very 
high. This indicates that for these sources possibly a 
two or multi-temperature greybody spctrum is more 
appropriate. 
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Observations at 1300 /im by Mezger et al. (1988) 
have revealed 6 protostellar condensations within an 
area of 3' x 3' around NGC 2024. The far-infrared 
maps presented here being of limited angular resolu- 
tion do not resolve these condensations. However the 
205 yum and also the 100 /im maps show evidence of 
north-south extension with the northern tip having a 
slight westward bend (position of FIR-1, 2 ; Mezger 
et al. 1988). Positionally B9 coincides with FIR- 
5 (brightest condensation detected by Mezger et al. 
1988) and also with LBS 33SM, a CS core detected 
by Lada, Bally, & Stark (1991a) (hereafter LBS), later 
detected at 1300 /im by Launhardt et al. (1996). 
The HIRES map at 100 /im shows possible signs of 
anomaly, particularly near the peak. The flux density 
measured (in 3' dia) from the 100 /im map shows a 
very low value of ^18 kJy compared to both the IRAS 
PSC value (35.3 kJy) and results from the KAO obser- 
vations by Thronson et al. (1984) (85 kJy). We have 
used the measurements of Thronson et al. (1984) in 
the Table 3 as well as for determining Td and (3 for B9. 

Figure 6 shows the dust temperature distribution 
for this region. It has been generated using the same 
method as for Orion A. It shows several isolated peaks 
with temperatures varying between 60 to 70 K for most 
of them. The source B3 is quite cold with a maximum 
of 30 K and going below 15 K towards the edges. The 
T150 map (Figure 6) shows very interesting structure 
close to the main source B9. It distinctly shows that the 
region around B9 has an elongated structure with indi- 
cations of at least 2 condensations. This structure is 
very similar to the structure shown by the NH3 column 
density map generated by Schulz et al. (1991). This 
very well substantiates the fact that these are indeed 
protostellar condensations since they are so luminous 
in the infrared as well. The cold source B3 is found to 
be mostly optically thin (peak 0.07) at 150 /im. 

The total luminosity of the region calculated us- 
ing the same method as for Orion A (Section 3.1) is 
1.0 xlO 5 Lq. The diffuse emission contributes about 
50% (at both the wavelengths) of the total luminos- 
ity. This is substantially larger than the value measured 
in Orion A (Sec 3.1) and in other young Galactic star 
forming regions (~ 35% in W3 1 ; Ghosh et al. 1989) . 
The morphology of the [CII] (158 /im) emitting region 
in Orion B (south) as observed by Jaffe et al. (1994) 
agrees very well with the maps presented in this paper. 
The total [CII] luminosity detected for their assumed 
distance of 415 pc is 23OL , while we have estimated 
the total far-infrared luminosity for a distance of 450 



pc. Correcting for this difference in distances assumed 
we find the ratio of the [CII] luminosity to the total FIR 
luminosity to be 3 x 10" 3 . 

The observations presented here, being on a large 
scale and of good angular resolution would be useful 
as guidelines in understanding the star formation sce- 
nario vis-a-vis the temperature distribution for future 
large scale sub-millimeter continuum observations. 

3.3.2. Region around NGC 2023 

Figure 7 shows intensity maps of the region sur- 
rounding NGC 2023 at 205 and 138 /im alongwith 
HIRES maps at 100 and 60 /im. The maps at 205 and 
138 /im are generated by masking NGC 2024 prior to 
MEM deconvolution of the larger map so that details 
of the fainter source NGC 2023 are better visible and 
comparison with 100 and 60 /im HIRES maps is easier. 
The source NGC 2023 appears to have a double struc- 
ture at 205 /im. It is of interest to note that the map 
of fluorescent H2 observed by Gatley et al. (1987) ex- 
hibits a shell structure with main peaks separated by 
similar extent as B13 and B14. The intensity map at 
138 /im however shows only a single source close to 
the brighter source (B14) at 205 /im. Diffuse emission 
extending southward from the main source is clearly 
visible at 60, 100 and 205 /im, while at 138 /im the 
limited dynamic range of the map restricts this emis- 
sion to a seemingly isolated peak. Table 3 presents the 
positions, flux densities and associated sources for the 
2 peaks (B13 and B14) detected in this region. The 
coordinates of the global peaks in the 60 and 100 /im 
maps is found to be - 5 h 39 m 7 "- 2° 17' 44". These 
global peaks are separated from B14 by more than 1 '.5. 
Since there is no other detected source of such bright- 
ness in the neighbourhood and observations at 100 and 
138 /im trace dust of similar temperature we have as- 
sociated B14 with NGC 2023. Both sources B13 and 
B14 have been detected at 1300 /im by Launhardt et 
al. (1996), but there is no IRAS detection of B13. 
For both sources we have used the available flux den- 
sities at wavelengths longer than 60 /im to estimate Td 
and (3 from a single temperature greybody spectrum 
(Table 3). The flatness of (3 for B14 indicates the ex- 
istence of dust components of at least two temperature 
ranges. The source B13 is found to be reasonably cold 
(Td =16 K), thus justifying it's non-detection in the 
IRAS wavebands. 

Figure 8 shows the T( 138/205) and T150 maps for 
the region around NGC 2023 (derived assuming (3 = 
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2). The T( 138/205) map shows a peak which does 
not coincide with the intensity peaks at 205 and 138 
fim. The T( 138/205) peak, however coincides with the 
peak positions at 60 and 100 /Ltm. The temperature of 
most of the region is around 30 K and in places goes 
down to values as low as <~ 20 K. The T150 map more 
or less traces the high intensity regions. The T150 dis- 
played for regions which lie outside the contours of 
T(138/205) map are for an assumed temperature of 10 
K (not plotted in the T( 138/205) map). 

NGC 2023, being a well known reflection nebula 
with a photodissociation region (excited by HD 37903) 
has been studied in various molecular and fine struc- 
ture ([CII] at 158 11m) lines. The continuum maps pre- 
sented in Figure 7 are morphologically similar to the 
C ls O map by Wyrowski et al. (1997) and CO (J=2- 
-1) map by White et al. (1990). The [CII] observa- 
tions at 158 fim by Howe et al. (1991) and Jaffe et al. 
(1994) detect emission which is approximately 25 % 
of the emission from NGC 2024. The rise in temper- 
ature, T( 138/205), to the north-west of the main peak 
is explained by the presence of the exciting star HD 
37903. The total luminosity measured in this region is 
1.5 xlO 3 L Q . The diffuse emission contributes about 
35% of the total emission, which is smaller than the 
contribution of diffuse emission around NGC 2024. 
The total [CII] luminosity from the region as detected 
by Howe et al. (1991) is 7 L Q and is equal to - 0.5% 
of the total FIR luminosity estimated here. This is con- 
sistent with results from most of the Galactic star form- 
ing regions as given by Howe et al. (1991). Detection 
of low temperature in the region around NGC 2023 
alongwith emission in the far-infrared, support ongo- 
ing formation of low to intermediate stars in the region. 
Indication of more than 50 % of the emission being in 
the condensations further indicates star formation in 
clumps in this region. 

4. Summary 

In this paper we have presented large scale high an- 
gular resolution (<~1') far-infrared (205 and 138 fim) 
maps of the Orion A and Orion B (south) molecular 
clouds. For limited regions, maps at 57 /xm have also 
been presented. HIRES processed IRAS images at 100 
and 60 /im of the Orion B have been used for com- 
parison. The regions covered include OMC-1, 2, 3 in 
Orion A and NGC 2024 and 2023 in Orion B (south). 
In all 15 condensations in Orion A and 14 condensa- 
tions in Orion B (south) have been identified. In OMC- 



3 the very cold source MMS 6 (Chini et al. 1997) 
was detected at both 205 and 138 fim and the tempera- 
ture was determined to be ~15±2 K. The far-infrared 
counterparts of the extremely dense CS cores (Lada, 
Bally, & Stark 1991a) LBS 33, LBS 34, LBS 36 and 
LBS 40 were detected in Orion B (south). The to- 
tal far-infrared luminosity and the relative contribution 
of diffuse emission to the luminosity of these regions 
have also been quantified. 

We thank B. Das for his help during initial stages 
of analysis of the 57 [im data and S.L. D'Costa, M.V. 
Naik, D.M. Patkar, M.B. Naik, S.A. Chalke, GS. 
Meshram & C.B. Bakalkar for their support for the 
experiment. The members of TIFR Balloon Facility 
(Balloon group and Control & Instrumentation group), 
Hyderabad, are thanked for their roles in conducting 
the balloon flights. Thanks are due to IPAC for provid- 
ing HIRES processed IRAS data. We thank the anony- 
mous referee for the suggestions which have improved 
the quality of the paper. 
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A. Estimation of dust temperature distribution using the two band FIR observations 

As mentioned in Section 2.1 the spatial intensity distribution in each band is gridded with pixels of size 0.'3x0f3. 
Only those pixels which have signals exceeding 5 times the noise level of the corresponding map are considered to be 
"valid" for the purpose of finding the dust temperature distribution. Each intensity map is first smoothed by taking a 
running average over 3x3 "valid" pixels. This is done to be extremely conservative about the derived spatial structures 
in the temperature maps. Since the two intensity maps refer to identical portions of the sky and the beams are also 
almost identical, a "pixel by pixel" ratio between the two maps gives reliable dust temperatures. 

For optically thin emission the observed flux density at any wavelength can be written as 

F v = QB v (T d )T v (Al) 
where £1 is the solid angle of the region under consideration. 

We assume r v oc . Since the pixel sizes are identical the ratio of flux densities at any two wavelengths would be: 

Fv^ = fll\ 3+0 [exp(hv 2 lkT d )-l] 
F V1 \V2j [exp{hv\/kT d )-\] 

Thus for an assumed emissivity exponent (/3) the ratio of the flux densities at any two wavelengths is a function 
of the dust temperature only. The temperature is computed using a lookup table containing the calculated ratios of 
flux densities for various temperatures and an assumed /?. This process is repeated for all pixels to generate the dust 
temperature distribution map. 
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Figure Captions 



Contour levels are same as in Fig. 7c. 



Fig. 1. — Intensity map of Orion B (South) at 205 fjm. 
Dotted boundary shows the region mapped. The peak 
is 2.50 kJy/sq. arcmin. The contour levels are at 90 
to 30% (in steps of 20), and 20, 10, 5, 2.5, 1.25 and 
0.63% of the peak. 

Fig. 2. — The intensity maps of Orion A. Peaks are 
9.00kJy/sq. arcmin at 205 /mi and 29.1 kJy/sq. arcmin 
at 138 /im. The contour levels are same as in Figure 1, 
with an additional contour at 0.31% of the peak. The 
lowest contour level is 3 times the measured noise. 

Fig. 3. — The intensity map of the Orion A at 57 /mi. 
The peak is 59.6 kJy/sq. arcmin The contour levels are 
same as in Figure 2. 

Fig. 4. — The dust temperature map for the region 
covering OMC-1 and 2, assuming a dust emissivity of 
e v oc v 2 . The lowest temperature found in this region is 
14 K. The numbers written on the contours refer to the 
temperatures, with the highest contour displayed being 
at 70 K. 

Fig. 5. — The intensity maps of the region around 
NGC 2024. Peaks are (a) 2.50 kJy/sq. arcmin at 
205/mi, (b) 7.53 kJy/sq. arcmin at 138 /mi, (c) 3.80 
kJy/sq. arcmin at 100 fjm and (d) 19.1 kJy/sq. arcmin 
at 57 /mi The contour levels (a) and (b) are same as in 
Figure 1. Levels in (c) are 60 to 10 % (in steps of 10) 
and 5, 2.5 and 1% of its peak. Levels in (d) are 90, 75, 
60, 40, 20, 10, 7.5, 5, 3.5, 2, 1 and 0.7% of the peak. 

Fig. 6. — (a) Same as in Fig 4, except being for the 
region around NGC 2024. (b) T150 map for the same 
region. Contour levels are 90 to 10 % (in steps of 10) 
and 5, 2.5 and 1% of its peak (0.24). 

Fig. 7. — The intensity maps for the region around 
NGC 2023. The peaks are (a) 325 Jy/sq. arcmin at 
205 /mi, (b) 329 Jy/sq. arcmin at 138 /xm, (c) 511 
Jy/sq. arcmin at 100 /mi and (d) 786 Jy/sq. arcmin at 
60 /im. The contour levels in (a) are 95%, 90 to 10% 
(in steps of 10) and 5% of the peak and in (b) are 95%, 
90 to 30 % (in steps of 10) of the peak. The contour 
levels in (c) and (d) are the same as in Fig. 6b except 
for an additional contour at 95 % of the peak. 

Fig. 8. — (a) T(138/205) map for the region around 
NGC 2023. The highest contour displayed is at 40 K. 
(b) t 150 map for NGC 2023 with peak value of 0.14. 
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Table 1 

Position, flux densities (in 3' dia), associations, T d & f3 of sources identified in Orion A at 205 

AND 138 



Source # 


Position 2 (1950) 
R.A. Dec. 


Association b 


^138 

kJy 


^205 

kJy 


^1300 

Jy 


K 


p 


Al 


5 31 51.6 


- 5 06 22 


IRAS 05318-0506 


1.26 


0.58 




29 




A2 


5 32 02.0 


- 4 56 30 




1.24+ 


0.96 




20 




A3 


5 32 04.6 


-5 19 14 




0.83 


0.39 




28 




A4 


5 32 11.6 


- 5 02 28 




0.56 


0.36 




22 




A5 


5 32 32.8 


- 5 22 57 




2.65 


0.68 




80 




A6 


5 32 43.7 


- 4 57 42 


IRAS 05327-0457 


2.12 


0.96 




30 




A7 


5 32 47.3 


-4 58 33 




2.06 


0.86+ 




>32 




A8 


5 32 47.7 


-5 24 51 


OMC-1 


112.50 


35.70 




68 d 


0.5 d 


A9 


5 32 55.6 


- 5 04 00 


OMC-3,MMS 4,5,6,VLA 3 


0.42 


0.53 


4.77 


15 


2.0 


A10 


5 32 59.5 


- 5 14 25 


FIR 6a,6b,6c,6d,VLA 14 


2.22 


1.44+ 


7.45 


25 


1.6 


All 


5 33 01.4 


-5 11 16 


FIR 2,3,4, VLA 11,12 


1.88 


1.50+ 


6.94 


20 


1.9 


A12 


5 33 01.6 


-5 13 01 


FIR 4,5,6a,6b 


1.98+ 


1.71 


10.4 


20 


1.8 


A13 


5 33 07.1 


- 5 01 49 




0.43+ 


0.54 




15 





+Not detected as a local peak in this waveband. Flux density is obtained by integrating in a circle of 3' dia 
around the coordinates tabulated. 

a Units of right ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 

b IRAS : IRAS PSC ; FIR & MMS : Chini et al. (1997); VLA : Reipurth, Rodriguez, & Chini (1999). 

°Single temperature greybody spectrum; if flux densities at only 2 bands are available emissivity e„ oc v 1 
is assumed. 

d Using F 57 , Fi 38 & F 20 5 (this work) and F 450 = 19.8 kJy, F 790 = 2.671 kJy & Fn 00 = 0.989 kJy (Goldsmith, 
Bergin, & Lis 1997). 



Table 2 

Position and flux densities (in 3' dia) of sources identified in Orion A at 57 /xm 



Source # 


Position 3 (1950) 


F57 


Associations 




RA. Dec. 


kJy 




1 


5 32 47.3 - 5 24 27 


219 


OMC-1 


2 


5 33 00.2 -5 16 18 


5.91 


M43 


3 


5 33 03.3 -5 21 46 


3.99 





a Units of right ascension are hours, minutes and seconds, 
and units of declination are degrees, arcminutes, and arc- 
seconds. 
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Table 3 

Position, flux densities (in 3' dia), associations, T d & (3 of sources identified in Orion B at 205 

AND 138 LlM 



Source # 


Position 3 (1950) 


Association b 


F12 


F 2 5 


F57 


F60 


F100 


Fl38 


F205 


F1300 








R.A. 


Dec. 




Jy 


Jy 


kJy 


kJy 


kJy 


kJy 


kJy 


Jy 






NGC 2024 




























region 




























Bl 


5 38 21.4 


- 1 50 19 


05383-0150 


60 


65 


NM d 


0.95 


1.40 


0.67 


0.29 




32 


2.5 


B2 


5 38 21.9 


- 1 54 24 




68 


68 


NM 


1.17 




1.04 


0.30 




25 


3.8 


B3 


5 38 43.0 


- 1 49 55 


05387-0149 
LBS 40 SM 


80 


133 








1.69 


1.13 


2.6 


20 


2.3 


B4 


5 38 44.5 


- 1 58 44 




87 


120 


4.50+ 




2.79 


2.18 


0.67 




21 


3.8 


B5 


5 38 52.6 


- 1 46 10 


05388-0147 


89 




NM 






0.57+ 


0.49 




20 




B6 


5 38 57.2 


- 1 54 00 




110 




5.51 






1.89 


0.90 




95 


0.4 


B7 


5 38 59.4 


- 1 47 52 




81 


175 


NM 


2.83 


3.59 


1.25 


0.37+ 




27 


3.9 


B8 


5 39 08.4 


- 1 52 13 


05391-0152 


248 


1060 


6.62 


10.74 




4.24 


1.31 




32 


2.8 


B9 


5 39 13.7 


- 1 56 52 


05393-0156 
LBS 33 SM 


86 


916 


71.70 




85.00 e 


32.30 


12.06 


115 


59 


1.0 


BIO 


5 39 21.8 


- 1 45 47 


05392-0145 


17 


31 


NM 


0.32 




0.83 


0.59 




35 


0.5 


Bll 


5 39 21.8 


- 1 49 13 


05393-0150 


80 


176 


0.71 + 


2.10 


2.96 


1.56 


0.65 




32 


2.4 


B12 


5 39 28.5 


-2 01 36 












0.77 


0.51 


0.27+ 




80 


0.2 


NGC 2023 




























region 




























B13 


5 39 02.9 


- 2 20 16 


LBS 36 SM 


25 


52 


NM 






0.79+ 


0.82 


4.7 


16 


2.1 


B14 


5 39 14.6 


- 2 19 08 


LBS 34 SM 


14 




NM 






1.93 


1.41 


12.5 


25 


1.3 



+Not detected as a local peak in this waveband. Flux density is obtained by integrating in a circle of 3' dia around the coordinates tabulated. 

a Units of right ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes, and arcseconds. 

b LBS SM : CS cores (Lada, Bally, & Stark 1991a) also detected at 1300 /im by Launhardt et al. (1996). IRAS PSC for all other entries. 

c Single temperature greybody spectrum; if fluxes at only 2 bands are available emissivity e v oc v 2 is assumed. 

d Not mapped at 57 /im. 

e Thronson et al. (1984). 
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' " Orion B (South) "\ : ..r -; 



. I I 1 1 L 

5 h 41 m s 39 m 40 s 38 m 20 s 

Right Ascension (1950) 




Right Ascension (1950) 



Right Ascension (1950) 




5 h 33 m 40 s 33 m s 32 m 20 : 



Right Ascension (1950) 



-5° 6' 




5 h 33 m 30 s 33 m 15 s 33 m s 32 m 45 s 32 m 30 s 



Right Ascension (1950) 




Right Ascension (1950) 



Right Ascension (1950) 




Right Ascension (1950) Right Ascension (1950) Right Ascension (1950) Right Ascension (1950) 




Right Ascension ( 1 950) Right Ascension ( 1 950) 



